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ABSTRACT 


Thxs thesis makes a study of »' Yagi-Uda 
Antenna”. The design of an Yagi array, given the 
gain, array size and percentage bandwididi is given, 
A study of current distribution in the elements of 
the Yagi, by Moments Method and also by Eing's 
method is explained. The procedure for calculation 
of current distribution, gam, input impedance and 
directivity by Ring’s method are given, a computer 
program was prepared, which analyzes the Yagi by 
King’s method. Experimental set ups were used to 
measure the performance. The computed results are 
compared with earlier published results and with 

i 

experimental measurements. 
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INTRODUCTION 

1 1 General 

Antennas are basic components of any electrical 
communication system which depends on free space as a 
propagating medium. The antenna is the connecting lihR 
between free space and the transmitter or receiver. In 
many systems for navigational or direction-finding purposes, 
the operational characteristics of the system are designed 
around the directive properties of the antenna. In other 
systems, the antenna may be used simply to radiate 
energy in an omnidirectional pattern in ordei to provide 
a Broadcast type of coverage In still other systems, the 
antenna may have a highly directional pattern to achieve 
increased gain and reduced interferenoe, 

A general broadcast receiver antenna must have the 
following properties Low cost, simplicity of construction 
and erection, low wind rosistanoe, low sensitivity to small 
shifts inposition and must require least maintenance, 

This IS because the users are mostly non te clinical people 
Television chamois fall in tho very High Froquonoies 
(VHP) and Ultra High Preq.uenoi es (UHi') bands. Since VHV 
and UHP signals decay rapidly v,ith distance, TV receiving 
antennas are required to have a high gam, particularly 
infringeareas , 
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Yagi-^da antenna fxts the till adequately in all 
respects. Its structure has low wind resistance and its 
heamwidth is sufficiently large to take care of small 
shifts in position, due to wind, rain, birds etc It has 
a gain in the range of 8 to 1 3 dB over a simple dipole and 
the construction does not involve any costly machining 
process* Yagi's arc used in Radio astronomy applications 
also 

A conventional Yagi an ay consists of a row of 
parallel straight cylindrical dipoles of which only the 
second one is driven by a source and all others are para- 
sitic# The driven element is approximately half a wave 
length long and is usually the second in number The first 
element is slightly longer than half a wavelength and is 
known as reflector The remaining elements are known as 
directors which are shorter tlx-xn the driven element 
Biguro1*i1 shows tho parameters of a Yagi antenna 

The properties of an antenna which are most often 
of interest are tho radiation pattern, gam and impedance. 

In this thesis these aspects of Yagi antenna are studied 

1 • 2 Outline of the work 

This chapter gives an inbroduction to the suboect 
and presents a literature survey on Yagi antennas Chapter 2 
explains the design of Yagi antenna given some specification 
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FIG 11 ATYPICAL YAGI - UDA - ARRAY 
1 REFLECTOR 
? DRIVEN ELEMENT 
3- N DIRECTORS 


4 


like Directivity, Bandwidth, array size etc 
Tables given in the Appendix I aid the design. Chapter 5 
presents an analysis of this antenna by Moments method. 
Computer techniques more than analytical discussions is 
the important difference between this method and the 
King’s three term current theory described in tho next 
chaptei (4) King’s method, unlike his predecessor's, 
follows an analytical approach to the foimulatlon of the 
current expression. His matrices are of a smaller order 
and henco subsequent computer work is comparatively simple 

Chapter 5 03!plains the computations involved in 
finding the characteristics of the antenna like Input 
Impedance, Gain, Directivity and lield Pattern, Chapter 6 
gives the actual experimental set up for the above measure- 
ments Since tho author wanted to verify the claims made 
by Cheng and Chen (17) about increased directivity with 
opt imijaation of elements, the antenna built, conforms to 
their design specifications Chapter 7 presents a discussion, 
on the results obtained A computer -J.ioting appears in the 
Appendix I H which will aid any furth) r work on this antenna 

1 • 3 l iterature Survey 

Prom a recent paper by Kahn (1), it is confirmed now 
that Prof Uda of Tohoku Imperial University, Japan invented 
the Yagi-Uda array (2) in the year 1926 However as Yagi, 
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a collegue of Prof Uda, published a paper in English (3) 
in 1928, in which he first mentioned the array, the 
antenna was known as 'Yagi airay' Wow it is called 
rightly as 'Yagi-Uda array* 

In 1946 Walkinshaw (4) computed and plotted a large 
number of directive polar diagrams for several shortcffhd- 
fire arrays with a driven half wave dipole and a maximum 
of four half wave parasitic elements Reid (5) and 
Pishenden and Wiblin (6) analysed the array but their assu- 
mption that equal currents flow in all directors was wrongs 
Reid considered the array as a continuous arrangement of 
very short dipoles Rishenden and Wiblin considered only 
half wave elements with sinusoidal current distribution. 

Ehrenspeck and Poehler (7) investigated experimentally 
and systematically a method for obtaining maximum gain 
from a Yagi-Uda array with equally spaced directors of 
equal length They concluded that the phase velocity 
of the surface wave travelling along the row of directors 
could be used as a design criterion The phase velocity 
depends on the element length and spacing parameters in a 
complicated way The surface wave concept has also been 
used by various authors to calculate the phase velocity of 
infinitely long uniform dipole arrays with assumed current 
distributions (8), (9) to determine the cut off frequencies 
(10) and to analyze infinite and semi-infinite Yagi-Uda 
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arrays (11)» (12) An approximate method for studyxng the 
tehaviour of finite arrays hy matching the terminal Zone 
solutions of semi-infinite arrays was discussed hy 
Mai Hour (12) 

To maximize gain the antenna parameters e g inter- 
element spacing, length of the elements and diameter of 
the elements must he optimized This approach is found in 
Bojsen et al (13), Morris (14), Tseng and Oheng (15) and 
Cheng and Chen (18), (19) papers Boo sen's results dispute 
the travelling wave theory But he neglected the dipole 
radius and assumed sinusoidal current distributions* Moms 
used an arrangement in which the driven element was a half 
wave dipole ( 2 h 2 =v\/2) and the single reflector (21i^=0 51 A ) 
was placed at a quarter wavelength behind the driven 
element (b ^2 = ^ 25/s) This gave a high forward gain He 
plotted curves showing the variation of gam with inter- 
element spacing. His conclusion that if the interelement 
spacing IS longer than (0.4/x), the gam will fall concur® 
with the findings of Ehrenspeck and Poehler (7) 

Oheng and Chen used a three term approximation for 
current, developed by King et al (16) They maximized 
gam by optimizing interelement spacing (18) and element 
length (19) Kao fez found that nonlinear optimization 
reduces the side lobes of Yagi antenna (22) 
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Shen gave a n-umerical approach (21) for the design 
of Yagi antenna (20) Thiele (24) and Takla (23) also 
have contributed for the analysis of Yagi array* 

Thiele (27) and Harrington (17) have explained the 
moments method solution for current and field problems 
of Yagi array. 
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CHAPTER 2 

DESIGN OP YAGI-UDA ANTENNA 


2 1 General 

In this chapter the theoietical constraints on the 
parameters of the Yagi-Uda antenna are given* Also the 
design of an Yagi array, given a set of specifications 
like Directivity, Bandwidth etc is indicated The 
Tables given in the Appendix I aid the design 

2 2 Theoretical Constraints 

The parameters of the Yagi antenna are the number 
of elements^the length of the elements, the ladius of the 
elements and the spacing between the elements. 

In King's (16) 3 term current theory, the assumptions 
Ip^ob:^^ and are made to satisfy some conditions 

Here [3 q(= 2Tr/^ ) is the free space wave number and A is the 
wavelength, 'b' is the inter-elcment spacing and is the 
half element length of the ith conductor in the array Also 
Po where 'a' is the radius of tho elements, 

Pishenden and Wiblin ( 6) stipulated that tho length 
of the array should not exceed SAand 'b' should be A/3 
Poehler (7) Morris (14) and Cheng and Chen (18) found both 
theoretically and esjperimen tally bhat the spacing between 
the reflector and the driven element, if kept eq.ual to 
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0 25 A gives a gain i^lso the Reflector should 

be 0 51 A long Inter director spacing must be less 
than 0,4 A 

Gheng and Ohen found that in practice if the length 
of tho driven element is slightly less than 0 5A , the 
gain IS increased 

The Bandwidth of an array is defined as the frequency 
range in which the directivity falls to 5 dB less than, 
the maximiM directivity value. According to Shen (20) the 
bandwidth of an Yagi array is limited to the bandwidth 
of the pass bands of the travelling wave Yagi array is 
supposed to support a travelling wave since it radiates 
primarily in the end fire direction The pass bands are 
bounded by the upper and lower cut off frequencies (^Qh is 
essentially a frequency term) The upper cut off frequency 
can be found by a rigorous mathematical analysis of tho 
structure (28) The lower limit is determined by the 
requirement that the power density should decay at least 
as fhst as exp (~r/h) where 'i' is the distance measured 
from the centre of the array to a point in the transverse 
direction. The three pass bands are 

Kh = 0 to 1.57 
Kh = 4.49 to 4 71 V 
and Kh « 7.73 to 7 85 ! 


Kb 3.14 
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2 3 DesifSyi of an Yagi array 

Sh. 0 n (7) has indicated a procedure by which, given 
the central frequency, required Bandwidth and the length 
of the array, we can design an Yagi antenna which would 
have maximum directivity and minimim number of elements. 

The design constraints wore that the central frequency 
should bo 200 MHz (/b= 1.5 m) and the bandwidth approximately 
10% The length of the array should be 5 metres 

The length of the array being 2Aand tho Bandwidth 

10%, we chose curve 'b' where = 1 0 for our design in 

Figure 2 1 This gives a directivity of approximately 

in Appendix I 

1 2 dB. Now referring to Table ( II )i which pertains to ^ - 

b/h = 1 , we can see that, if the number of parasitic elements 
IS equal to 9, we get a directivity of approximately 12 dB 
with 10,3% bandwidth The jf^^h value 1 36 determines 
the length of the directors approximately as 0 65 m As 
b/h ^ 1, the spacing between tho elements = 0 324 m ¥e 
can take a/h =0.01 and then a = 0 0032 m The array length 
will be 1 96/\or 2. 9^1 m. Tho array will have one reflector, 
one driven element and eight directors. 

2,4 Conclusion 

With the help of the graph and the tables wo can 
with 

design any array even/^only two out of the throe constmints 
namely Bandwidth, Directivity or array size given. 


CHAPTER 3 


analysis or yagi antema by mombfts method 

3,1 Gen eral 

High speed digital computers have simplified the 
3 oh of the antenna engineers in that numerical techniq.ues 
can he used where analytical solutioreare not avaxlahle. 

It IS now possible with the aid of approximation techni- 
ques to get very accurate results about antenna performance 
with the aid of computers. 

Field computations hy Moments Method was first intro- 
duced hy Harrington (17) Thiele (27) solved the Yagi 
antenna problem using Moments method In his work Matrix 
methods are used largely in solutions to field problems. 

The basic idea is to reduce a functional equation to a 
matrix equation by known techniques 

5 2 Method of Moments 

Let us see a general procedure for solving lin^r 
equations, called the method of moments. Consider the 
inhomogeneous equation 

L (f) = g (5.1) 

where ’L' is an operator, 'g* is the source or exitation 
(known function) and 'f* is the field or response 
(unknown function to be determined) We are interested 
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xn an unique solution Let 'f be expanded in a series 

of functions f^,f 2 ,f 2 *• in the domain of L as 



where theo{j^ are constants. 1/e shall call the f^^ expansion 
functions or basis functions. Por exact solutions (32) 

IS an infinite summation and f^^ form a complete set of 
basis functions. For approximate solutions, however, (32) 
IS a finite summation. Substibuting (3 2) in (3.1) and 
using the linearity of L, we have 

f "n (fj = S (3.3) 

It IS assiMed that a sin. table inner product 
has been determined for the problem Defining a set of 
weighting functions or testing functions, w^ , W 2 » • 

in the range of L, and taking the inner product of (3.3) 
with each w^^^ we get 

"n n ^ '^n^ ~ (3.4) 

m = 1,2,3, . This set of equations can be written in 

matrix form as 


L = C 


where 



L f.j^ J' 1 ^2^ 

L f.j^ <W2> L ± 2 / 



(3 5) 

(3 6) 
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, g> 
^W2, g> 

t 

• 


(3*7) 


If the matrix f IJ is nonsingular its inverse ^ exists* 
Ihecx^ are then given by 


= LCI r^' (3.6) 

and the solution for f is given by (3.2) Defining 

the matrix of functions, 


l^rJ = L^l- ^2- *3 ■ 1 (3 9) 

we write, 

W= = 1 «nl i Cl LC (3.10) 

This solution may be exact or approximate depending upon 
the choice of the f^^ and Wj^ The choice w^^ = f^^ is known 
as Galerkin's method 


One of the main tasks in any particular problem 
IS the choice of the f^ and w^^* The f^ should be linearly 
independent and chosen so that some super position (3»2) 
can approximate f reasonably well. The Wj^ should also be 
linearly independent and chosen so that the product 
<(^Wn» g^ depend on relatively independent properties of g. 
Additional factors which affect the choice of f^ and w^ 


are 
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1 . the accuracy of solutions desired 

2, the ease of evaluation of the matrix 
elements 

3 the size of the matni that can be inverted 

4* the realisation of a well conditioned 
matrix C 1 J 

The integration involved in evaluating the 
^mn ” '^^m' ^n/* (5*^) often difficult to perfoim 

in problems of practical interest A simpler way to get 
approximate solutions is to require that equation (3*5) 
be satisfied at discrete points in the region of interest. 
This procedure is called point~ matching method 

There are two types of Bases fimction used in 
Moments Method ihe Entire domain Bases function is 
defined and non- zero over the entire domain. Examples 
are Pouiier, Ohebyshev, Maclaurin and Legendre series 
The other types are sub-domain bases, which exist in the 
domain of Bqp? But are zero over part of that domain. 
Examples are triangle functions, quadratic interpolmtiojjs 
etc 

5 3 Applj cation to the Yagi-Uda Antenna 

Thiele has analysed Yagi-Bda antenna using Moments 
method. He used a current generator as a source on the 
driven element Entire domain Basis function and Point 
matching techniques were used. 
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Pxg 3 1 shows the coordinate system used to 
analyze Yagi-Uda array To find in 

%(/) = (?) 

where J is the current density and ^ is the incident 
electric field ¥o consider the various elements in the 
array as individual wire antennas ^ It can be shown that 
the operator for the pth element is 

f-^p/2 


IS a good choice (27) Here 


J G(r,r')dZ' (3 12) dz 



S (r,r'5 = (i+jkr)(2r^-3a^)+k^a^rf3 (3 13 ) 

and r J ( X' -X) ^+( Y» ~y )2 + (Z'-Z)^4a2 (3 14 ) 

Since we are considering each element in the array 
to be an individual wire antenna let us use the entire 
domain basis functions 


The current Jp(Z) on the pth element is assiamed 


to be uniform along the circumference Thus 
Jp(Z) = 2F a Ip(Z) 


(3.15) 


We expand I ( Z) in terms of entire domain basis function 
as follows 


H “L L 

I (Z) I cos (2n-1) 

P n=1 % ^P 


(3.16) 




I 
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FIG 3 



COORDINATE SYSTEM USED TO ANALYZE 
YAGI UDA ARRAY 
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The operator L^^CJ) must be evaluated not only 
when the observation point is on element p, but also 
when it IS on the other elements as well, as indicated 
in Pig* 3*2* Thus we can write 


% 


' P=1 




h/2 

^ <J(r>r') 

-^p/2 


CO b( 2n“1 ) 


TTz' 


dZ' 


(3 17) 


where (D+2) represents the total number of elements in 
the Yagi-Uda array (D=number of directors ) and N is the 
number of entire domain basis functions retained on each 
element in the array. The entire field radiated by the 
array is represented by (3*17), 


In an array composed of a reflector a driven 
element and D directors, let us assume that there will 
be N modes on each element and let each element be of 
different length Using (3*1 7), the first part of the 
system of equations is then of the form 


D+2 







^m ,n 

= 0, 

m = 1,2, .NID (3 18) 

P=1 

n=:1 

* p 

P 


whe re 

it 

Ai^/c 

S/2 

$ 

-^p/2 

G-(r,r') cos(2n-1) dZ’ 

P 


(5 19) 
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FIG 3 2 DIAGRAM SHOWING DISTANCE FROM A MATCHING POINT 
ON pth ELEMENT TO SOURCE REGION ON qth ELEMENT 
INSERT SHOWS RELATIONSHIP BETWEEN z' AND 0 WHEN 
OBSERVATION POINT AND SOURCE REGION ARE ON THE 
SAME ELEMENT 
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These equations are generated by requiring that 
tangential B field be zero at N points on each director 
Or E IS zero at NxD points on the directors The 
matching points on a director are shown in Pig 3 3 

The next N equations are similar to the previous 
BxD equations since tangential E vanishes at N points on 
the reflector element, as shown m Hg 3 3 Thus 
D-i-2 N 

.ad-. ^m.n ^n ~ m = (Nxl))+1, Nx(D+ 1) 
p=1 n=1 P P . ^ 

(3^20) 

The last IT equations are generated by using the 
boundary condition on the driven element as shown m 
Pig 3*4. That IS 
D+2 it 

.§1 -SZ In = 0, m = H1(D+1)+1, 

P=1 I* P m(Ii+2)-1 (3.21) 

and for the driven element at the driving point 

21 = 1 , e = D+2, m = IT i^(D+2) (3.22) 

n =1 ® 

On the driven element the tangential E bounda 3 :y 
condition is only enforced at points even though 
there are N modes The Nth equation on the exiter 
arises from the constraint on the terminal current value. 

In equation (3*11) we have sot (iP") to aero as no 
incident field is considered in this formulation Let us 


z 


N 

N 1 

li 

\ 


flG^i PARASITIC ELEMENT WITH 
I OINl S Al ONG ITS AXIS 


amp^ 


N 1 
■N 2 


Z4 

Zi 

^2 


FIG 3 4 DRIVEN ELEMENT WITH N 1 
POINTS ALONG ITS AXIS 


N - MATCHING 


MAI CHING 
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consider an example given by Thiele in which one 
reflector, one driven element and one director are 
present Assuming that only two modes are present m 


each element 

, we 

get 
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That IS at point on element 1 (the reflector), 
^ IS the ''field" generated by mode 1, and Z^ g 
field generated at the same point by mode 2 is the 

field generated by mode 1 at point Z2 on the reflector 
and Z22 IS that of mode 2 Similarly we can explain 

» ^ 54 * ^45 ^ 44 * 


If we write ( 5 * 25 ) in a submatrix form, 




* r 





S12 

^13 j 1 

h'. 

) 

^ 2 . 


“0 

j ®21 

^22 

^23 ' 

> 


0 


^32 

* 

33 i 

I3 1 


0 



—1 

j 

I ^ 


1 




4 


j 






0 ' 

1 

1 






0 , 


( 3 . 24 ) 
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then, it IS clear that regardless of the niamler of sut- 
matrices, those on the main diagonal of suhnatrices 
( Si 1 ^ 22 * ^33^ will represent uhe field generated by the 

current on the element at which the tangential B boundary 
condition is being enforced Hence if there are D 
identical directors the last D submatrices on the main 
diagonal will all be identical 

Foi elements off the main diagonal of submatrices 
the explanation is as follows, The quantity represents 
the field at a point 1 on element 2 (the driven element) 
due to the first mode (mode 5) on the director (element 5) 
Similarly Z^g, and Z^g, Thus submatrices 1 

represent the interaction between elements q and r. If 
all directors arc of the same length then, = S_„ for 
director subnatrices Burthor lor uniform director 
spacing there will be several int erdirector spacings or 
distances that will be the same II the directors are 

also of same length, then wo need not calculate 

the director submatrices individually Thus computation 

time IS saved. 

Once the current in the individual element are 
found by solving (3*24), the calculation of pattern, 
gain and impedance are very simple The procedure is 
explained in Chapter 5* 
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3.4 Conclusxon 

As Olieng and Chen (Id) pointed out, in sub- 
sectioning the array eloaents, matrices of much larger 
dimension would have to be manipulated Thiele, 
himself said that two modes per element is very inadcq.ui,te 
and minimiJm three per parasitic elements and five 
modes for the driven element are needed for good results 
This limits the nimiber of airay elements that can be 
handled Furthermore, the currents on the parasitic 
elements depend much more critically upon those on other 
mutually coupled elements The effects of small errors 
multiply and theio would be convergence problems unless 
more subsections than those normally required for 
driven elements are taken. 

In using the three teim theory of King's, the 
largest matrices to be handled for an N clement array are 

JF 

of a dimension IMf and no convergence problems are 
encountered in optimizing tho array for maximum gain 
Henco in this work. King's three term theory is used 
for current distribution and the procedure ig explained 
in the next chapter. 
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ANALYSIS OP YAGI ANTMNA BY KING'S METHOD 


King, Mach and Sandler (l6) have formulated 
an alternative approach to the problem of finding 
current distribution in the elements of the Yagi-Uda 
antenna. They assumed a three term approximation 
for current which is given below (for the ith. 
element) 


I,(ZP=AA (4.1) 

1 \ Zj. 


z 


= 0 for X fk 2t A 2 > and D^ are complex constants 
to be evaluated and 


^o = Sin P^Ch^ - \Z^\) ( 42 ) 

^1 

Pq = cos - cos p^h (4.3) 

^1 

Hq = cos \ p^Z^ - cos ^ p^h^ (4.4) 


These three components have been given a 

physical explanation by King as follows M is the 

°Zi 

simple sinusoid This component of the current is 
maintained directly by the generator ^ the driving 
voltage m the driven element number 2), It does not 
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include the components that are induced by coupling 
between different parts of the antenna The currents 
induced by the interaction between charges moving 
in the more or less widely separated sections of the 
antenna appear in two parts. One of these the shifted 
cosine, IS maintained by that part of the interaction 
that is eQ.uivalent to a constant field acting in phase 
at all points along the antenna The other part the shifte 
shifted cosine with half angle arguments is the correction 
that takes account of the phase lag introduced by the 
retarded instead of instantaneous interaction, 

let us try to determine the integral equation 
for base driven monopole (Tig. 4«l) over perfectly 
conducting ground screen (which may be approximated to 
a cylindrical dipole) We have the boundary conditions 
E 2 (Z) =5 0, on the surface F = a of the perfectly 
conducting antenna. 

Prom Maxwell Lorentz 8q.uations we know that 
for time dependent fields 


V X 

B = 

-V/ 

i‘o ® ) 

(4*5) 

V 

B = 

a/ 

0 

(4.6) 

y X 

B = 

3sjoB 

(4.7) 


/V 



r . 

B = 

V 

= 0 



ll 


6-Generntor 
1/2 V« 6(2) 



FIG 4 1 (a)COAXI AL LINE TERMINATED IN OPEN '^ND 

(blBASE DRIVEN MONOPOLE OVER PERFECTLY 
CONDUCIING GROUND SCREEN 
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whieh govern the interaction of charges and currents 
on conductors in space k convenient method of solving 
the vector partial differential equations (4.5‘-4*8) 
is with the use of scalar and vector potentials 0 and A 
It IS well known that the fields are derived from the 
vector potentials by the following equations 

B = V X A (4 9) 

and B = - V0 - j u)k. (4 10) 

f\J 'V-' 

If we substitute (4 9) and (4*10) in the Maxwell's 
Equations we get mixed vector equations for A and 0 . 

r 

The variables can be separated if the following conditions 
relating A and 0 are imposed. 

V A = --'f 0 (4.11) 

This IS known as he lorentz condition The 
resulting vector Helmholtz equations for A and 0 in 

nJ 

air are 

(xy2 ^02)^ ^ j ^2) 

+ Pq^) 0 = = 0 (4 13) 

As we consider thin cylindrical conductors all aligned 
in the Z direction, we see that the boundary condition 
B^(z) =s 0 on the surface a, leads to the following 
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B = - V 9 jA gives V 9= 

.c'. 

and "V A = ■“"* 9 gives 9 = "''-r ^ ^ f^om wliich 

P (■ v/ 


o 


■99 = -^o y ( ^ A ) Since has only a z-component 

p «- j 

on the surface and this is a function of z alche, we get 


d i. 


2t 


d z 


+ p 


"z 


0 


d z 


(4 14 ) 


on the surface P = a 


Equation (4.14) has the general solutions 


\{Z) = (0^ cos p^Z+a 2 sin p^Z) 

(4 15) 

where the symmetry conditions (-Z) = 12 ( 2 ), AgC-Z) = 
A(Z) are imposed, 0^^ and <^2 arbitrary constants of 
integration, c is the velocity of light. 


If we use froe space Green's function to 
derive the particular int^rals of (4.12) we get for 
A 2 for a thin cylindrical conductor of length 2h and 
radius a with its centre at Z = 0 as 


A. 




.0. 

411: 


A 


h -3^0® 

1^(2') 2-5 — az' 


-h 


(4.16) 
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/ 2 2 

where - free space wavelength, and R=Y(Z-Z ' 

For a point on the surfaceU ~ a, we combine (4.15) and 
(4.16) to get the integral equation for the current as 


, h -3^0® 

= i I^cz-) az' 

° -h 

= ^^[Oioos PqZ+Oj sin 1 Z{ ) 

^ 0 

(4.18) 

where = Vu-/e_ = 120ic ohms, 

»•> V vj \J 


The one dimensional lorentz condition is 


"2) Z “ 

C. 

(4.19) 

From (4 15) and (4.19) 

kn 

_5_ ) 

0 a 


H (z) = ( 

Po 


= -CJj^ Sin PqZ+ ^ 

+ d^cos p Z for Z > 0 
- CJcos p®Z for Z < 0 



(4.20) 


Etsr definition the driving voltage of 


the delta function generator is 

0 [0 (z) - 0 (-Z)] = v^® 

Hence 0^ == | 


(4,21) 
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4jL Agf ( 2 ) — C 
-h 


-3PoR 

l^iZ') ^ g dZ' 


^ [o^Qos PqZ+ 2 sin \ 2i ] 

( 4 . 22 ) 

In the above equation (4.22), the constant 0^ oan be 
removed by the following manipulation , consider the 
following expression 

4Tt [ A2(Z)~A2(h)] = K^(Z,Z’) dZ ‘ 

-h 

= C^l 2 pQ ^ ^ ^ 

(4 25) 

which is obtained from (4.22) by subtracting At^ A^C^i) 
from both sides. 



-h 


and K^(Z,Z’) = K (Z, Z » )“^(h,Z’ ) (4.25) 

The constant 0^ oan now be expressed interms of 
U and by setting Z=h Since the left side of (4.25) 
then vanishes the right side can be solved for 0^ to give 

-(I 4- u) 

^ COS P^h 


(4 26) 
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If this value of 0^ is substituted in to 
( 4 , 25 ) we get 

h D 4 t^ T 

f \{i') VZ.Z-) dz' =^~rp;h [a C Po(h-lzi) 

+ TJ(cos PqZ - cos P^h) ] 

(4*27) 

This IS an exact expression and does not involve 
any approximation. 



The current l 2 ;(-i^) vanishes at the ends Z = jyh 
and is continuous through the generator at Z=0^ and 
satisfies the symmetry condition I^(-Z) = I^(Z), The 
Kemal in (4»22) is 

-oPoR f 

K (Z, Z’) = ® , R = f(Z-Z')^+a^ (4-28) 

Separating into real and imaginary parts 

cosp.R sinp R 

Kj^(Z,Z‘)= ,Kj(Z,Z') = 

( 4 . 29 ) 

The dimensionless quantities E^(Z,Z*)/Pq and 
Kj(Z,Z’)/Pq are shown graphically in Rig, 4.2 as 
functions of p^\^Z-Z'/ Ej^(Z,Z’)/Pq has a sharp high 
peak at ZssZ’, Its magnitude (l/p^a) is very large 
compared with 1 since if has been assumed that 1 




0 12 3 


po( Z~Z ) 

sin OoR cosOoR 
__ — , — and cos 


FIG 4 2 THE FUNCTIONS 
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On the other hand (Z,Z')/p^ varies only slowly with 
p^jZ - Z‘) and never exceeds the value 1 Also (sin p^E/p^R) 
IS very well approximated by cos (p^R/2) in the range 
0 rS Pq 1 2 - The value of cos (p^R/2) is hardly 

affected if the small quantity a is neglected and 
P^R IS approximated by Pq\ Z - Z(. 

Consider now the real and imaginary parts of the 

3 P ^ 

integral ) I^CZ') --- dZ’ = P (say) 

“h 

P = + jPj (4.30) 

h *^o^ 

Then Pj^(h,Z) = 5^ 12,(2') dZ' 

= (z)i(z)=y^ i(z) (4 31) 

P^ (h,Z) = - C l2(Z') dZ' 

-h 

= -. l2(Z'’)cos I p^(Z-.Z')dZ’ 

-h 

(4.32) 

The reasoning for (4.31) is as follows. Since 
the Kernal is quite small except at and very near 
Z' " Z. where it rises to a very large value, it is 
clear that the current near Z'=Z is pranarily significant 
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in determining the value of the integral at Z, Or 
simply, the integral is approximately proportional to 
l(Z), The proportionality constant is best 
determined where ^2^^^ is a maximum 

This can he further modified by considering that 
the integral on the left of (4.31) becomes luite small 
at the ends of the antenna where Z = ± h, the RHS 
vanishes identically at these points since (±h) = 0, 
Hence a better approximation of (4.31) is 


h 


PR(Z)~Pj^(h) = r l 2 (Z')[J 2 f^(Z,Z')-K^(h,Z ')3 dZ' 

-h 


= % 

where '|^2 a new constant^ 


(4.33) 


Equation (4,32) can bee further simplified 
h 


Pj( h,Z)= ^ I 2 (Z‘) cos ^ p^(Z-Z»)dZ’ 


•h 


(4.32) 


h 


--Po 1 l2(Z’)[cos ^ P^(Z^Z^)f-3^ 


cos I p^(Z+Z’)]dZ» (4.34) 


h 


2 Pq cos ^ p^Z ( l2(Z')cos i p^Z' dZ’ 


(4.35) 
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For antennas that do not greatly exceed 
P^h = Tt xn electrical half length, specifically for 
p^h ^ 511/4, 

Bj(h,Z) = Pj (h,0) cos ^ p^2' (4*36) 

where P^(h,0) = -2 Pq ^ 2 ^o^' 

o 

4*4 Dnlvation of Three ~ Term Fo r mula 

ITow in equation (4*27) if we introduce the 
approximations (4 33) and (4 36), and consider the 
i.H S. of (4*27) 

h h 

ms ^ { l 2 (Z»)K(Z,Z')dZ‘ - ^ I^CZ')^ (h,Z’)dZ' 

-h ~h 

(4.37) 

Since K^(Z,Z') = K(Z,Z») - K(h,Z’) 

h 

-a2 p^( J i2(z>) 00s I P^Z' dZ')ooe 1 P^Z] 
0 

-[0-d2Pq ( ^ 12 ( 2 ’) cos ^ p^Z’dZ’) cos 1 p^h] 
o 

(4,38) 

= y 2 Po^°®® 2 ^o^) 

o 


(4 39) 
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Also U = — ^ ( I 2 (2') IC(h,Z«) dZ' 


■h 


( 4 . 40 ) 


* 3 ^ h. 

[0-02 P„ 00 s i p^h I^CZ') 

O 

COS ^ Pq Z« dZ»] 

( 4 . 41 ) 

■4;*' [3 1 P„ 00a i p^h I2 (Z') 


00 s I p^z' az'] (4.42) 

Oonsiderxng the RHS of ( 4 . 27 ) 

™s=/^;h[2 V Po " 

/ O’^o - 1 ^ 

(cos P^Z-COS P^h) (32 P^ cos JP^h) C 12(2’) 


cos ^ P^Z' dZ'3 (4 43) 

Now transferring the second term of IHS (4,59) 
to RHS (4.43) and equating IHS and RHS we get, 

'-'0 '■^O 

h - 

+ [ r I„ (Z') COS ^ P Z'dZ'].[Tp~“ . 

0 , ' 

3‘^2p 

cos I p^h+ 23 p^(cos ^ p^Z 


(cos p^Z-cos p„h) . 

'^0 0 ' /\.Tt 


cos ^ P^h)] 


(4 44) 
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This can be written as follows 


12 ( 2 ) = A[sin p^(h-| Z\ ) ]+B[ cos p^Z-cosp^h] 

+ D[cos 2 Pq 2 - cos 2 


where 


0 2ii V 

A ^ 


e 


O^oosp^hfj 


(4.45) 
(4 46) 


B 


'T 2 cos ? H^h h 

J a_2_ (2.) 

1 2 cos u h ■> 


and 


D = 


2 00s ^ 


cos * P^Z» dZ’ 

2 ^ O 

a2P, h 


(4.47) 


'V 2 ^ 2 


I^,(Z»)cos ^p^Z» dZ’ (4.48) 


Thus we have got a three term approximation 
for current. 


The above formulation holds good for the 

Single thin oylindrioal conductor case This theory 

can be extended to two conductors (see Appendix II) 

i 

and by induction the general case of N thin cylindrical 

/V 

conductors arranged along the Z axis, lagi antonna is 
a special case of this array in which N-1 conductors 
are parasites and only one element is deriven. The 
solution of this array we will see in the next section. 
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4 • 4 Application of Three term Theory t o Yag i ar ray 


On the hasis of the three term approximation 
the current in the single driven element has the 
form 


(Zg) - ^^oZ2'^^2 ^oZ 2 ^2^oZ2 
The currents in the parastic elements are 

= Vo2l + ^ = 1-5.4 » 

(4 4913 ) 

where 

“o2K Po(‘"K - 

^oZK =°°® h ~ °°® Po % ^4 50) 

HoZK = 2 Po^K - =°® 2 l-'o ^IC 

and are complex coefficients 

and mufst he evaluated in terms of the 
driving voltage The integral equation for the 

driven element is 



( 4 . 51 ) 
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Here 


_ JS _ jS- 


R 


Kx 


R 


Ka 


(4 52) 


where 


^Ki ~ + ^li» %ih= /(VZ,’)^+b| 


and 


% = -4r^ 


1' -ICi 

(4 53) 


h. 


1=1 “h. 


) <52-' 


(4.54) 


and ^ Q = 120u: 


The integral equations for the remaining (H-I) 

parastio elements are 
iig 

^ j ^oZ'2 ^42 d ^^2 

-h2 

N 

N 

5- 21 \ Hozx' ’%d ^hc'K > ^2’ 

1=1 


a4iL 


^cosp^hg % ^oZIC .N 


(4 55) 
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In order to obtain approrimate solutions of the 
N simultaneous integral eq.uations, use may be made of 
the properties of the real and imaginary parts of the 
iCernal That is 


®oZ'k: bccds ' ''oZ’c (4 56 ) 


where stands tor ^oZ'lv ^oZ'K 

^^iQCdR ^ real part of the Kernal, 


^oZ'k (4.57) 


For any distribution Henoe 


A ]^k 


»'liOT = ( M (Zp Z', ) azj 


7itav **ozE ■•■ '■t’Lcav ®oze 


(5 58) 


Xlk 

=5 ^oZ'iv (2 t^,ZjP dZj 


'xQaiJ ^oZK ■^^lacdu ^ozic 


(4.59) 
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Irt 


-h. 

= V 


I 

£ 


IffidD ^oZiC +'<^WdD 


h 


(4,60) 


where 


stands for ’’equal to by definition” 


-f* 

The term fiQjjij, F^zj; is added in order to provide 
greater flexibility and symmetry %en i k and p^b ^ 1, 

It has been shown^-ftp^S^^^^C) 

h-i 

-h *^^1 ^oZlC (4.61) 

1 

^1 

^oZ’i ^Sidl (2j^»2^) (4.62) 

Therefore 


"JCIT - y P^2£ + fkiav HoZJf 

(4 63) 

= l^idir ^oZk +vldU \ZK (4.64) 

= H-T^iD r„2E: +'1' h^^j. C 4 .S 5 ) 

so we can write from (4.58 - 4 SO) and (4 65-4 65 ) 
and (4. 51) 
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m ii 

^ [V 22dV ^oZ2 ^oZ2 ^ 

N f h. 


\ 2idU ^oZ2 ^^■’%LdU ^oZ23 

1=1 

N f h 

+ ^ t'^2idD ^oZ2 f^aU-dD ^oZ2^ 

1=1 jk. 

= f ^SFTh: t 2 ^O^OZS-^ U^0Z2] 

o o 2 


( 4 . 66 ) 


Prom (4.55), the (N-l) equations are 
■^t^^I^2dV ^oZK "^^'cPdV ^oZIv^ 


Y 


h 


ICidU ^^oZlJ 


N f 

+ :SL. ®i[^KidU ^oZK 
1=1 

N f h 

+ ^ ^f'^KidD ^oZK ■^'^KidD ^oZlJ 
1=1 


“ cosp h.. "oZX 

- O ^0 ii 


Uir P. 


x{=l,5,4. ^ 
(4.67) 


These equations will be satisfied if the 
coefficient of each of the three distribution functions 
IS individually required to vanish i.e in (4 66) 



(4.69) 
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hi V 22dV y 2idy 2idD ^ ~ ^ 


(4.70) 


Similarly in (4.67) 




- ^ Uir° 

1 0 


(4 71) 


^ ^2dV ^ f®i^ I0LdU‘''\ ^KidD^"° 


(4.72) 


(4 73) 


Defining ICronecker 6 as 

6 _ J 0 1 k 

^iZ " 11 1 = K 

we can combine (4.69# 4 70, 4.71 and 4.72) 
The 2N equations are 


h Z2dV t^^iV CidU ^ °i ZidD ] 

1=1 


3 47111 

Po^K - ~er~ 

~y 0 


- 0 


(4.74) 


■^VKPdV t^WzidU kidD^ " ^ 

1=1 

with K = 1,2,3 Z, These equations together with 
(4.68) determine the (21T+1) constants h* 

1 = 1,2,3, N 

Now to evaluate the function Ujr, we define 
J "" ^oZi» Zi^(V^i ^ 


(4.75) 





4-5$ 


XX 

i " I'oz'i (+-V6) 


-h 

1 

■where 

r^^o^Eah 


(4 77) 




Kih 


Siah = ^ 


? 

IxC 


From (4.54) 

,r -J 
^K = 471 


R 


1=1 


Ki V ^ ^XV^ r^u ^ ^^i"^ 


iCiD 


(4^78) 


Since only element 2 is driven, A^=0, 1?^ 2, so 


I 

Urr = -Tr“ ^ V iQ(jy (hj^)+ ^ [B^ r 


ZL 

'K ^ 471 ' f 


1=1 

\ "^IClD ^ (4 79) 

The substitution of (4 ^7) xn (4 71) gives 


C Vl^2av *^o^iC 3 

N . f 

^ ^ 'iviU^^iv^" 't'jtidU 

N 

\ ■'H^KidD 

(4 80) 



With 1,2, N let 
0J2V == 

(4 81) 

^I&.u= -V^,iaa =0® Po*^iC <'^•82^ 

^IilD “ ~ l''EidD “°® ^oV (4.8J) 

With this notation (4 >)t ) with (4,74) gives the 
following set of 2N equations for determining the 2E 
coeffioients and B in tp>rms of 

^ \ ^KiD ~ ^I2>V ^2 

1=1 

^ t KidU \ '^V^dD ~ ^dV h 

(4.85) 

with E = 1,2, h 

Here 0 |j, 0 ^^, and \^i aie(HxH) matrices and 

02y» ^ 26 . 7 * ^ ^ (Hxl) vectors. 

The matrix form of ( 84 ) and ( 85 ) are 


[^)tj]|Bj+ [ 0^] (4.86) 

[Yaff] + [4 an] '■ »'(= -h 
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4.5 Evaluation of 4'*- 


Since each integral in (4 58-4 60)^ (4 63-4 65) 
and (4 75-' 4 77 ) is approximated by a linear conbination 
of two terms with arbitrary coefficients, these can 
be evaluated by equating both sides in ( 4 . 58 - 4 . 60 ) and 
(4 63-4.65) at two values of Z The values chosen are 
Z=0 and 2 =hj^y 2 addition to Z=h^^ where both sides 
must vanish, can write 

= V it -54 


“h. 


M 


oZ 


, 3 _ (o,zp az^ (4.88) 


¥ (bs) 


A, 


( l.r (Z*) K, .(-f f Z') 
) Vi' 1' lud 2 ’ 1' 

-h 


dZ, 


h 

( 

-h 


f ''oz'd %a(-2 > 4 ) ^4 ('t 


h 


=" \ I ^UiC2|) %,(0 ,Z') dZ^' 


-h 


Ul^"l^ "^ICLd' 


h 


i"- ^'oz-d %af°’4> 

-h. 


(4 90 ) 




= i " ^oz-x laa 91) 


-h 
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Wt 




-h. 


h 


= ( H^z'x \xa 
h 


(4,92) 


''iCxD^^’ = ''x’ 5 ''^Dx ) ^ICld 2“’ “i 


-h. 


f" ^'x> “1 


Here K = 1,2,3 .H, These integrations can be converted 

into generalised sine and cosine integral functions 
Then a high speed computer can calculate the ’ ¥'s. 

Once the W’s in (4# 88- 4 93) have been obtained for all 
values of i and K, the coefficient V may be determined 
from the eq,uations (4 58 - 6l) and (4*63 - 4»65) A.t 
2=0 these becomes 


‘f'TacdV ^ luidv 

(4 94) 


(4.95) 

^iidU^^”°°® ^o^K^'^'*'lidtrt^“°°®^ °2 

(4 96) 


(4.97) 



4 *? 


kt Z - 




“^JCKV^r^^ 


( 4 . 98 ) 


VlidV 1^ ^o^K ^■^H^Kidvl^°°^ “4 


B h. h. 

- cos )J= ^ 1 i4 K 


(4 99 ) 


P ^il-ir h. .r 

- cos -^^— ]= ^Ciu^ 2^ ^4 100) 


f h. ^ o^k 

'^KidD [°°® PoV 2 -°°® Po*‘El+‘t-'>CiaB[ = °® 4 - 

00 s -|-i ]= ) (4.101) 

Th.e solution of these eq.uations for the V are 
obtained directly They are 

^ KICdV = ^ L cos ] 

~ ^ 1-cos j (4.102) 

(4.103) 
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V lo.d'V ^2 2^^^ [1-cos - ^T^jiO) 

[cos cos 1 K 

(4 104) 

y KiDV 2 ^ - oos ] 

- ] j 1 1 (4 105) 

y iCidU ~ ■^2 — OjiA j 

" %U^^C/2^ ij (4 106) 

h _ -If' / \ 

V iCidU -4^2 {^(10^2 ^ [l~Qos - ^ICiU^C)) 

=0® PoV]} (4 107) 

V JiLdD “ -*2 [°°® - oos-'^^— ] 

■ [l-=os ]j (4 108) 

^Ha.aD“ “^2 [ l-=os p^h^^.] - ¥jj^^(0) 

r ^n^lr -r 

[ocs cos p^h^]? ( 4 .iog) 

Where 

Ai = sin p^h^^[cos ((^o^i{:/4^"°o®(Po^i(/2^3 

sin 1 ^^”®°®^*^o^i{:/2 5I (4.110) 



and 




A>2 = [1-cos p^hj^][G 0 S -. cos 3 

P P o^Z 

-[cos "^2 "■ Po^F^f 1-cos — ^ ] 

(4 111) 

Thus first all the ^ are determined from 
(4 102 to 4*111) Then the 0 matrix elements are 
computed from (4*81 - 4 83)* Finally solving 
equation (4 'b^^) 

1 1^] [ «Sd1 
Faul ^di)] 





- 027'^2 
1“ f 2dV^^2 


112) 





¥e get B^’s and D^'s was already calculated from 

(4.68), Substituting all these in ( 4 , 49 ) and (4 50) 
gives current in all the segments of Yagi-Uda airay. 


Am* No. A 63J1J 



CHAPTER 5 

COHPUTATIONb POR PATTBBH, &AIH AND IMPEDANCE 
5 1 Geneial 

Once the current in each segment of every 
element of the airay is found, we can calculate the 
Input Impedance and Admittance of the Yagi antenna. 
Radiation field pattern^ gain of the antenna and 
directivity of the antenna using computers 

5*2 Input Impe da no e 

The input impedance of the Yagi antenna is 
measured at the feed points of the driven element i»e 
at the center of the 2nd element in the array The 
input impedance helps us in designing the receiver 
circuits 

The input impedance calculation in a computer 

analysis is very simple As we assume that the 

exitdtion IS one volt, the reciprocal of the current 

at the middle of the second element, I„ (21,2) will 

^2 

directly give us the Input admittance The reciprocal 
of the Input admittance i>o the input impedance 

5 3 Pattern Measurement 

The far field or radiauion field of an antenna 
is one of its most important characteristics The field 
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paxtern is actually a three dimensional or space 
pattern, and ito complete deociiotion requires field 
in I ensity calculations in all directions in space* 

A space pattern can he calculated according 
to uhe following procedure Let us jceep the antenna 
at ohe origin with the x plane horizontal and the 

u 

Z a^is vertical as in Pig 5*1 -i-’hen on an imaginary 
sphere of large radius with the originat the centre, 
patterns of 0 and 0 components of the electric field 
(1 )q, 00) are calculated along latitude circles (that 
IS circles of constant latitude or polar angle, 0)* 
These patterns are calculated as a function of the 
longitude or azimuth angle 0 Calculating such 
patterns at 5° intervals in latitude from 9 - 0 to 
180° and 0 = 0°to 360^ completely descrihes the 
radiated field 


From King (16) we knox-r that the radiation field 
for a thin cylindrical conductor of length 2h and radius a 
with its centre at Z=0, 




3 Qh 


h 


Sin 9 r I^(Z’) — dZ’ 

"<-h 


(5 1) 


where R is the distance from an arbitrary point on the 
antenna to the field point and is given in oerms of r 
and Z' by the cosine law (I'lg 5 2) 
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+ (Z')" - 2 rZ* cos © (5 2) 

2 2 

In the radiation zone r (Z') If the binomial 
expansion is applied to (p 2) and only the linear term 
in Z' IS retained, uhe following approximate form 
IS obtained for R 


R = r - Z> cos ©, ( Por)‘^>> 1 (5.3) 

The phase variation of exp (~o p replaced 

with tne linear phase variation given by (5 3)i,o,by 
exp(“0 P 0 ^ ^ ^ o^' amplitude 1/R Of 

exp (-3 PqR)/R is a slowly varying function of Z' and is 
leplaced by l/r, where r is the distance to the centre 
of 1he antenna Since r is independent of Z' all 
funcbions of r may be removed from the integral in 5 1 


and the final form for is 




G 


•OPo’- 


sin © 




h 


IgCZ') c 


-h 


3 P Z' cos 9 


^5 4) 


dZ' 


Now to determine the electric field maintained at 
distant points by the currents in the N elements of the 
Yagi Uda array, we find thai for 

^ Z 2 ^^ 2 ^ ~ ■^2 Pq ^^^2 ~ 1 ^ 2 ^^ 

+ B 2 ( COB p Z^ - CO s p hg ) 

+ Dg (cos Z2-’00s ^ Pq hg) (5 5 ) 
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and 


(Z^) = (cos - cos Po \) 

+ (oos ^ Po \ - cos ^ Pq h^) 


the electromagnetic field is 


3 M* 

E (E2,S,0) = sin 


10 2 (5-6) 


w -3 P „ R. 

^ X e ° ^ 

1=1 Ri 


h 


iz ® 


3 pQ cos © 


dZ:^ 


(5 7) 


Referring to Pig 5*3, wo make the far-field approxi- 
mation R^ = R^-'fj^in amplitude However in phase, 

“ (1-2) h sin 0 cos 0 (5 8) 

where i is the ith element ntunher and 'h' is the 


\ - ^2 


inter-element spacing As we have defined d^= -(i-2)b, 
the equation (5*7) becomes 


(9.0) = 

r i-Q 


h_ 


1 0 Po'^i ^ ^ 

1 ® 


i h (q) 


e 


1=1 

3 P 0 cos 9 


sin 9 dZ' 


(5 


9) 
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Equation. (5 9) was used to calculate the 
E field pattern E('TT/2,0) c^xves the Azimutal or 
Horizontal pattern 

5 4 G- ain 

The gain of an antenna is defined as 

G-q-iyi _ Maximum Radiation i nt ensit y( te st antenna) 

~ Maximum Radiation intensity [Iseferenc e"" 

antenna ) 

(5 10 ) 

It IS ideal to consider a reference antenna which 
radiates uniformly in all directions with gain of 
unity But in practice as no isotropic radiator is 
availahle, we use a dipole as a reference which has 
a gain of 1 64 over the isotropic radiator 

We are more interested in Directivity since 
Yagi antenna is a directive antenna 


Directivity = 


M aximum radiaojon i ntensity 
Average radiation intensity 


(5 11) 


If we take that the antenna radiates a maximum 
in the direction (9 q»0q) ^ 

The radiation intensity U is related to the 
Pointing vector P which in turn is related to the 
Electric field intensity E(25) 




(5 12) 
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•where xs the intrinsic impedance of the free space 

Ihe ladiation intensity is the power per -unit solid 

2 

angle and is eq.ual to r times P Moreover by multi- 
plying the niimerator and denominator of (5 11) we get 
4 u |{ B (0 ,0 )/2 

Directivity = 

powex ladiated 


If ire take that the antenna ladiates a maximum in 
the direction (©q’^o^ to fcal power radiated 

2r 7^ 

= 1 d0 ) I E' (©,0)j ^ sin 9 d© d0 (5 14) 

0 6 

Hence Directivity is calculated from the radiation 
pattern and from (5 12)and (5 15 )t 

5 5 Oonclu sion 

The theoretical calculation of impedance, l field 
paotein and gain are described in ohis chapter A 
com-nuter listing appearing in the appendix will indicate 
the calculation procedure 


OHiLPTER 6 


ANT®TEA MmSUREMBNTS 


6.1 Ssjissal 

In the last chapter we saw the theoretical 
calculation procedures using a computer. In the 
present chapter the practical measurement set up is 
explained 

6.2 Mtoaa 

As we designed the Yagi antenna for channel IV 
whose centre frequency we took as 65 MHz, the following 
are the dimensions of the antenna for a initial set up. 
Since the wavelength (P\) is 4,615 meters, the radius 
of the elements is taken as 1 5 cms, the reflector 1 177 
meters, the driven element 1 150 meters and the 
directors 99 2 oms The spacing between the reflector 
and the driven element is 1 154 meters and the other 
interelement spacings are 1 45 meters 

However, in order to check Cheng and Chen's 
optimization results, the antenna was made with provisions 
to vary the element spaces as well as the lengths 
A locknut arrangement was used to move the eletneoits 
along the cross boom and concentric tubes were used 
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to vary the lengths of the elements The mam cross 
"boom length is 8 1 meters The height of the antenna 
IS approximately 3.5 meters from ground level 

k folded dipole was used for the driven 
element as this facilitates the use of the 3^0 O, tv 
cable as feeder* (Polded dipole impedances; 300 -n-), k 
balun was used to provide matching between the balanced 
and unbalanced circuits of the antenna The Pig 6*1 
shows the Baiun construction details (29) 

Por a test transmitter antenna, a folded dipole 
whose length was also 4.615 meters was constructed* The 
photoes at the end of the thesis explains the physical 
set up 

Experimental measurements on antennas become 
difficult on account of the fact that all antenna parameters 
are affected by reflections from nearby objects. To get 
over this difficulty, all measurements were made on the 
roof of the building in an environment relatively free 
from reflections, 

6 3 Imn edan oe, Measur.ement 

The impedance measurement set up is shown in 
Pig 6,2 A hewlet Packard Tector Voltmeter was used 
to measure directly the reflection coefficient (both 


i 


FIG 



T V Ribbon 
cable 

Balanced line Unbalanced line 

6 1 BAIUN CONSTRUCTION 4 TO 1 IMPEDENCE 
TRANSFORMER 






FIG 6 2 IMPEDENCE MEASUREMENT 
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67 


have to be off set by 30° in the opposite or OW 
direction, in order to correct for this initial phase 
off set. 

After replacing the short with the antenna, B/a 
1+1^ IB now measured to be 1,62 / +14° The point 
P must then be rotated 0¥ This angle can be plotted by 
striicing off from point D the distance SO with a 
pair of dividers to where it intersects the circumference 
of the chart at E, The magnitude of P of the antenna, 

OP IS then measured off from 0 along OE This point X 
equals 50 (5*5 + 0 1 0) or (275+0 50 ) ohms whose magnitude 
280 ohms is close to 300 ohms for a folded dipole, 

^•4 P attern Measurement 

To measure the far field pattern of the antenna, 
the distance between the Transmitting and the Receiving 
antennas must be appropriate. If this distance is 
too small, then the near field or Presnel pattern is 
obtained. The Presnel pattern is a function of the 
distance at which it is measured. Por accurate far-field 
measurements the antenna under test should be illuminated 
with a plane wave front Since plane wave fronts, 
are obtainable only at infinite distances, some limits 
must be specified* A commonly specified criterion is 
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that the phase difference between the centre and edges 
of the antenna under test shall b e no greater than A/16, 
If this IS the case, then from Fig, 6,4 

(R+6)^ = (6 1) 

and 

R^ + 26 ( 1 + 6^ = R^ + (6 2) 

2 

For R>>6-!^vd, 6 may be neglected, and 

a & (6,3) 

For 6 = A/16, R^ 2 d^/A (6.4) 

The receiving antenna aperture is approzimatoly 2 A and 
hence R A 8 A • We kept a distance of approximately 
10 ^ between the Transmitter and the Receiver (Fig, 6,5)* 
Moreover in order to avoid ground reflecti(3ns, the 
transmitter and receiver antennas were raoimted on two 
sections of the terrace of Western laboratory of 
IIT Kanpur with a valley between the two sections. 

The Yagi antenna is mounted on a Mast of 

height about 10 feet . At the base of the ma^t, 

an outer tube of a larger diameter is embedded in a 

square concrete structure, A circular disc with degrees 

marked on it for every 10° is attached to this outer 

has 

tube Tho inner tube (mast) /a pointer attached to it 



Transmitting antenna 
J aperature 

Rcc ei V inq antenna 
ap erature 


FIG 6 4 PHASE OlFFtRENCE Bf^TWtEN CENTER ANU 
Efrr 3F RECEIVING APR/W 
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so that when the mast and hence the antenna is rotated, 
the pointer indicates thi. angular movement accurately 
The arrangements are shown in Figs 6 fir, 6*6 and 6 7. 

The transmitter was adQusted for maximum output 
power. The antenna pointer was fixed at a reference 
angle 0 = 0°, The reading on the Jarrolds Field 
intensity iiieter was noted. 

Field intensity for every 10^ was noted and 
plotted in a circular graph, for an initial arrange ent 
of Oheng and Chen The final arrangemeit was then made 
for the elements and spacings of the array and the 
measurement and plotting was repeated Normalized 
values were used in the circular plot 

6,5 Gain Measurement 


The gam measurement was done by tho substitution 
method First a folded dipole was kept in the place of 
Yagi array and adoustments were made between the trans- 
mittei and receiver such that the receiver output is maxinium* 
Then tho dipole is replaced by the Yagi antenna. Again 
maximum Field strength is noted Then, 


Directivity = 


' R QqQiv..qy ■■ 
output of Receiver with Dipole 


(6.5) 




Lock nut 


FIG 6 6 FOLDED DiFOi E 




2 36 m 



1 - 6 Initial array patmq 
l' 6 Optimum array '^pacing 


FIG 6 8 PHYSICAL LAYOUT OF THE YAGl ANTENNA 
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The experiment was repeated with Kanpur T,Y, 
station as a transmitter Also gain for optimized 
array was measured The results of this chapter 
appears in the next chapter. 


1 


OHAPTER 7 


RESULTS km DISaUSSIdTS 


7.1 General 

There are three lypes of results which are 
available to us. King has published some results in 
his book ’Arrays of Oylmdrical dipole’ and Cheng and 
Ohen have given optimization results in their Technical 
Report The author of this thesis has formulated a 
computer program which appears in Appendix III which 
computes input impedance, current distribution in the 
elements, Directivity and gain and Electric field 
pattern The computer results are compared with King 
and Cheng and Chen results. The author has also 
constructed a Yagi antenna to the specifications of 
Cheng and Ohen to verify experimentally their claim 
of increased gam with optimum element spacing and 
length. The experimental and tneoretical results agree 
to a large extent, 

7.2 Computed resul t s vs Published Results 

King has given the current distribution in two 
full wave dipoles, one reflector and one driven element. 
Table 7 1 shows the correspondence between the published 
A 2 » and D^ values (the complex coefficients of the 


Table 7 1 


Item 

King’s Results 

Computed Results 


-0 249B - 04 

-OO 318E - 02 

-0,250E - 04 

-3 319E - 02 

®1 

0 708B - 04 

-dO 493E- 03 

0,752B— 04 

-3O 473B - 03 


0,ia3E - 03 
+0O.44IB-O3 

0 188E - 03 

+3O 437 B- 03 

h 

0.221E - 03 

+D 0, 456E'03 

0,217B - 03 

+oO,441B- 03 


0 439B - 03 
+0O.627S- 03 

0 445B - 03 
+3O 628B -03 


Table 7.2 


Item 


G-ain results by Oh. eng 
and Ohen(wr to isotr- 



as ratio as decible 


1.1 II I I I I ■ 11 II m nil — TataMU 

Computed 

Results 


1 

Initial 

array 

12,372 

10.92 

dB 

12,2, 10.87 dB 

2. 

Length 

perturbed 

16.42 

12.15 

dB 

16.08, 12,06 

dB 


array 






3 

Spacing 

perturbed 

19.16 

12,82 

dBB 

19 0, 12 79 

dB 


array 






4. 

Optimized 

array 

21 9 

13 4 

dB 

21 8 , 15.38 

dB 
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3 term currents distribution) and the computed values 
for a two element array 

Table 7 2 oompares the gain results published 
by Oh eng and Ohen with the computed values Their 
claim that optimization of length and spacing increases 
gain IS thus verified The corresponding impedance 
values also show that the program is reliable I'lg. 7*1 
shows the current distiibution in various elements of 
the array. The difference between the computed pattern 
and Cheng’s pattern published is negligible • 

7 3 E xperi me nta l Res ults vs Qompu t ed R e sults 

7 3 1 Eajdiaiiaii-Eatisiin 

Here the experimental set up described in 
Chapter 6 was used to see the optimization results. 

Fig 7*2 gives a polar graph representation of the 
pattern which also compares the experimental and 
theoretical values The experimental pattern shows lower 
side lobes and a higher back lobe than the theoretical 
curve This may be due to the fact that in theory 
the driven element is a single dipole but in the 
experiment, the driven element is a folded dipole Also 
the diameter of the elements were not exactly that of 
the theoretical value Still we can see that the peaks 
and troughs of the patterns coincide. 
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7,3 2 Input Impedance and Gram 

The impedance measurement values are compared 
in Table 7 3, G-ain values are also compared in the 
same table. Pig 7.3 oompareo the initial, opace 

perturbed and final pautoms. 

Apart from these experiments which were based 
on the published results, cesaputer analysis was carried 
on to find the Bandwidth of the antenna The computed 
values are shown in Table 7.4 and plotted in Pig. 7.4. 

It is worthwhile to evaluate the performance 
of a given Yagi antenna over several T? channels. This xs 
particularly true in situations like the lucknow TV 
broadcasting on Channel IV (centre frequency 65 MHz) and 
a relay station at Akrampur relaying, the same programme 
on Channel V (centre frequency 174 MHz ) and many 
times the same antenna is used for both. For this purpose, 
the gam of the Channel IV antenna for higher freq.uenoies 
was also computed and this is shown in Table 7,5, 

The pattern for Kanpur Television (Channel V) on 
a Channel IV antenna was found to have an almost 
non directional pattern, both experimentally and by 
computation. Still the poroximity of the high power TV 
station explains the prevailing practice of using 
Channel IV antenna for Channel V broadcast. 
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Table 7 3 



Imnedance 




Item 

Gheng and Ohens 

Exptai 

Gheng and 
Ohens 

Bxptal 

Initial 

array 

94.71+074.79 

100+070 

01 12 372 

11.45 

Final 

array 

10, 29+0 6 10 

10+0 3 

05 21,9 

20 08 


Table 7.4 

Fre<l 

MHz 

Gain 

wr to isotropic out 

Gain in dB 


56 

2,7'^ 

4.36 


57 

3.075 

4 88 


58 

3,736 

5 725 


59 

4.753 

6.77 


60 

6 23 

7.95 


61 

8.21 

9.14 


62 

10 57 

10 24 


63 

13 0 

11.14 


64 

15.1 

11.79 


65 

16 08 

12.06 


66 

10 318 

10.13 


67 

4 012 

6 03 
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k study of the increase in gam as the number 
of directors in a Yagi antenna is increased has also 
been done and the results are shown m Table 7 6 
and plotted m Pig 7.5 These results also agree 
with the generally held views like the increase in 
number of directors need not necessarily increase the 
gam A sort of saturation seems to occur with 
higher number of directors, 

7.4 Conc lus ion 

A computer program which is versatile is developed* 
Any future analysis can make use of this program. The 
antenna set up also can be used for future research. 

The effect of the diameter of the elements as an experi~ 
mental parameter can also be taken up. Optimization 
with cost, array size, impedance etc, can b e performed 
with little modification to the computer program. 
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Table 7 5 


( 

1 • 

i 

Gam 

■wr to isotropic out 

Gain in dB 

65 

12 52 

10,98 

80 

2 726 

4 35 

95 

2 721 

4 34 

110 

2.749 

4.393 

12 5 

2 76 

4.41 

140 

2.49 

3 96 

155 

2.25 

3.5 

170 

2 40 

3.8 

185 

2.98 

4 75 

200 

3.77 

5 77 


Table 7.6 


Number of 
elements 

Gain wr to isotropic 
radiator 

Grain in 
dB 

2 

3.54 

5.49 

3 

6,66 

8.23 

4 

9 97 

9.99 

5 

10.73 

10.30 

6 

13.34 

11.2 5 

7 

17 74 

12.49 

8 

17 81 

12.5 

9 

18 95 

12.77 

10 

23 83 

13.77 
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A PPMDIX I 

DESI&N mBLES FOR YAGI i^RAY 
Table I b/h = 05 ,-^= free space wavelengtli 


N 


DD 

(dB) 

Bandwidth 

1 

iLrray 
size ( ) 

6 

1 35 

7 6 

20 7 

0 65 

8 

1.34 

8 7 

20 2 

0 85 

10 

1.32 

9.6 

19.7 

1.05 

12 

1.31 

10 3 

18,3 

1.25 

14 

1.30 

11.0 

16 9 

1.45 

16 

1.29 

11 4 

15.5 

1.64 

18 

1.28 

11.9 

14.1 

1.83 

20 

1.28 

12 4 

13.3 

2.03 

24 

1 27 

13.3 

12,6 

2 42 

28 

1.26 

13.9 

11.1 

2.80 

32 

1.25 

14 5 

9 6 

3 18 

36 

1 25 

15.0 

8 9 

3.56 

40 

1.24 

15.5 

8 1 

3.94 

44 

1 24 

15.8 

8 1 

4.33 

48 

1 23 

16 0 

6 5 

4.70 

52 

1.23 

15 9 

6.5 

5.10 

5o 

1.22 

16.2 

4.9 

5.45 

60 

1.22 

16,0 

4.9 

5.83 
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Table 

II b/ii = 

1 0 

N 


D 

(dB) 

Bandwidth 

/° 

Array 
size ( ) 

3 

1 38 

6 8 

12.4 

0 66 

6 

1 37 

9 9 

11 0 

1 3 

8 

1,36 

11 3 

10 3 

1.73 

10 

1.36 

12 5 

10 3 

2 16 

12 

1.35 

13 3 

8 9 

2 58 

14 

1.34 

14 0 

7 5 

2 98 

16 

1 34 

14 7 

7 5 

3 41 

18 

1 34 

15 1 

6 8 

3 8? 

20 

1.33 

15 6 

6 0 

4 23 

22 

1 33 

15 9 

5.3 

4 64 

24 

1 32 

16 3 

4.6 

5 05 

26 

1 32 

16 5 

4. 6 

5 45 

28 

1.32 

16 6 

3.8 

5.85 

30 

1 32 

16 4 

3 8 

6.28 
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Table III b/h = 1,5 


IT 

^0^ 

D 

(dB) 

Bandwidth 

i 

Array- 
size ( ) 

3 

1.38 

7.9 

5.1 

4.91 

6 

1,37 

11 1 

4.4 

1.96 

8 

1 37 

12.7 

4 4 

2.62 

10 

1 37 

13 8 

4.0 

3.26 

12 

1 37 

14 8 

4 0 

3.91 

14 

1 37 

15 6 

4.0 

4.56 

16 

1 37 

16 1 

3 7 

5.20 

18 

1 36 

16 5 

3.3 

5.86 

20 

1.36 

17 0 

2.9 

6.50 


All these values are valid for the jPirst 


pass band P^h = 0 to 1 57) and for a/h = 0,01. 



APPENDIX II 


TWO ELEMENT ARRAY 


.on for two elements 


The integral equation (4.2?) for the current in 
a single isolated antenna is generalized to apply to 
the two Identical parallel and non-staggered elements 
shown in Eig, A 1 It is merely necessary to add 
to the vector potential on the surface of each element 
the oontrihutions by the current in the other element 
Hence, for element 1, the vector potential difference 

IS 

^ \ f 

-h 


^ 2 \o ^1 ^ 0 ^ 

^ o '^0 

- cos P^h) ] (A.l) 

Similarly, for element 2 

li / 

= ( [Iiz(2') K2^^(Z,Z*)+l2z(2’)K22^(Z,Z')]dZ 

-h 

= [ 2 '^20 Po(b.-\Z|)402(cosP„a-oos p^h)] 

(A,2) 
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1 2 



FIG A) TWO IDENTICAL PARALLEL 
ANTENNAS 
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In these expressions 

W2-2') = -55^-- 

= iC^^(Z,Z') - IC^3_(h,Z') a^3) 

“^’•^0^2 “^^0%2h 

XT / 7. 7. I 'I — ■£————— _ ,S—- 1. I— 

^ “ Ri 2 Rl2h 

= Kj_ 2 (Z,Z') - IC^2^h,Z*) (A,4) 

with 

= /(Z-Z')^+a^, R^^ = / (h-Z')^+a^ (A. 5) 

« ■ • 

R^2 = f(Z-Z»)^+b^, R^2h = V(h-Z')^ + (A. 6) 

B: 22 ^( 2 ,Z') and ^(Z,Z') are obtained from the above 
formulas when 1 is substituted for 2 and 2 for 1 in the 
subscripts, 

The two simultaneous integral equations (A,l) and (A*2) 
can be reduced to a single equation in two special oases. 


(a) The zero phase seq.uence when the two driving voltages 
are identical so that the two currents are the same, 

(b) the first phase sequence when the two deriving voltages 
and the resulting two currents are eOual in magnitude 
but 180^ out of phase. Specifically for the zero -r 
phase seq.uence> 
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""lo = \o = = l2z(Z)=l2<°^ (z) 

(A.7) 

So that the eq.uation (A,l) and (A. 2) become 

(Z') (Z,Z‘) dZ' 

-h 


r i 0.-JV, a /i, _ I 17. \ ^ -i. tt(®) 




[ ^ ain p^(h -iZ\) + u' 

(cos ~ cos (3^h )] (A. 8) 


where 


U 


and 


(0) 


(Z,Z») = ^ 




h 


^ K° (h,Z») dZ» 

-h 


47t 3 "-Z 




R 


11 


R, 


12 


(A. 9) 
(A.IO) 


(Z,Z') = (Z,Z’) - (h,ZO (A.ll) 


Similarly, for the first phase sequence, 

he = -ho = 

(A,12) 

so that the two equations again become alike and equal to 

i'*' (2-2’) dZ-= l^srh f 2 

-h ® o 

+ (cos p^Z - cos p^h)] (A.I3) 
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wliere 


(li,2’)d2’ 

-h 

(Z,Z’) = r 


Z^^^(Z,Z’) = (Z,Z’) - E^^^(h,2') 


(A. 14) 

(iU15) 

(A*16) 


A point to note is that the two phase se<3.uencec, differ 
only in the sign in E^°^ (Z,Z’) and E^^\z,Z'). 

If (A. 8) can he solved for the zero- phase- sequence 
current (Z) and (A*13) for the f irst- phase- se<luence 

current currents and maintained 

by the arbitrary voltages \o and ^2 0 obtained 

simply by superposition. This follows directly if 
and are so chosen that 

= 2 t^^io + = 2 C^IO - ^2o] 

In this case, 

^10 = 1^20 = 

so that 

^IZ^^^ ~ = I^^^^ (2r) + 

I22(Z) = Iz^°’ (Z) - (Z) (A.19) 
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^3^ope3:*~fcxes of the integrals t ^ 

The two integral eQ.uations (A«8) and (it. 13 ) for the 
pliase-seQ.U9xce currents are formally exactly like the e^luation 
(4.27) for the isolated antema. They differ only in the 
kernels of the integrals on the left and in the definitions 
(A» 9 ) and (ii.l 4 ) of the functions IT, Bach of these is 
now the algebraic sum of two terms that are identical 
except that the radius a appears in the first term, the 
distance h between the elements in t he second, term* In order 
to determine the effect of this difference on the current 
convenient to consider first the two extreme cases when the 
elements are very close together and when they are very 
far apart* 

The two elements may be considered close together when 
p^b”^! and b^ h. In this case, since b satisfies substantially 
the same conditions as a, the behaviour of the integrals 
that contain b corresponds closely to that of the integrals 
that contain a. ^lihen the antennas are so far apart 

(po4';>^l,b that + h^ - p^b)<?^' 1, the contribution 

to the difference kernels (Z,Z’) by the term 1^2 

- (e ° - (© ^ I2h ^ very small since R-,;, 

and \ 2 h slightly. In this case, the principal 

part of the interaction between the currents in the two 
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antennas is included in the function or and 

the integrals on the left in (A. 8) and (A*13) are onljr 
slightly different from the corresponding integral for 
the single antenna. The interaction between the currents 
in the two antennas is approximately as if each 
maintained along the other a vector potential that is 
uniform in amplitude and phase, Accordingly, the current 
induced in each elanent by the other is distributed in 
a first approximation as a shifted cosine. This conclusion 
follows directly from the fact that the component of 
current associated with the constant part of the vector 
potential along the surface of the isolated antenna is 
distributed in this manner, 

Wien, the separation of the two elemaits is such that 

f ^ 

P^b '^1 but not so great that f b + h differs 
negligibly from p^b, the vector potentials maintained by 
the currents on the one antenna at points along the surface 
of the other differ significantly f rom one another in 
phase due to retarded a ction. The induced currents should 
then have two components, the one distributed approximately 
as the shifted cosine with half-angle arguments, the other 
as the shifted cosine. 

In order to verify the correctness of these conslusions 
the difference integral 



has been evaluated for = -n: over a range of values of 
extending from 0,04 to 4*5. The real and imaginary 
parts are shown in Fig, A»2 together with the three 
trigonometric functions, sin (oos p^Z+l) and 

cos i 6 Z, to which the sine, shifted cosine and shifted 
cosine with half-angle arguments reduce when |3^h = it. 

For convenience in the graphical comparison - (cos j3^Z+l) 
and - cos ^ P^Z are shown. It is evident from Fig, A, 2 
that the real part of the difference integral approximates 
9in PqZ when PqZ< 1, 1+cosp^Z when p^b ^ 1. On the other 
hand, the imaginary part resembles the shifted cosine with 
half-angle arguments, in this case cos ^ stll 

Values of p^b, • 

As a consequence of these observations, the following 
approximate representation of the integrals In (A. 8) and 
(A, 13) is indicated: 
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